Spin Correlations and Finite-Size Effects in the One-dimensional Kondo Box 
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We analyze the Kondo effect of a magnetic impurity attached to an ultrasmall metallic wire using 
the density matrix renormalization group. The spatial spin correlation function and the impurity 
spectral density are computed for system sizes of up to L = 511 sites, covering the crossover from 
L < Ik to L > £k, with Ik the spin screening length. We establish a proportionality between 
the weight of the Kondo resonance and £k as function of L. This suggests a spectroscopic way of 
detecting the Kondo cloud. 

PACS numbers: 72.15.Qm, 73.63.-b, 73.63.Kv 73.63.Rt 



Scanning tunneling techniques have recently allowed 
to observe the Kondo effect of a magnetic atom in an 
ultrasmall metallic box possibly providing a direct 
probe of the long sought-after Kondo screening cloud. 
The Kondo effect is characterized by a narrow resonance 
of width ~ Tff, the Kondo temperature, at the Fermi 
energy £f 0- It is intimately related to the formation 
of a many-body singlet state, comprised of the impu- 
rity spin and a cloud of surrounding, spin-correlated elec- 
trons, the so-called Kondo spin screening cloud. Its spa- 
tial extent is vital for the coupling between neighboring 
Kondo impurities in a metal and, hence, is at the heart 
of spatial magnetic correlations and ordering transitions 
in Kondo and Anderson lattices and also in Hubbard or 
t — J systems, which exhibit local Kondo physics, as has 
been demonstrated by the dynamical mean field theory 
(DMFT) treatment of the problem However, while 
the spectral and thermodynamic features of Kondo im- 
purities have been well understood 01, the structure of 
the Kondo cloud has remained controversial for a long 
time. Researchers have been looking intensively for ways 
of observing the Kondo cloud. These include the Knight 
shift 1^ and recently theoretical investigations of the per- 
sistent current or the conductance (Q| in mesoscopic 
Kondo systems. For about 25 years it was generally be- 
lieved, and in the 1990s supported by scaling arguments 
0, that the Kondo cloud is characterized by a single 
length scale, = fivp/Tpc- It is the spin coherence 
length, i.e. the distance traveled by a scattered electron 
with Fermi velocity vp, until the impurity spin (whose 
lifetime is h/Tx) flips. Although can reach almost 
macroscopic values {S^k ~ lO^fc^^^ for Tk = 1 K, kp be- 
ing the Fermi wave number), it has never been observed 
in experiments. 

Only recently it was realized that another length scale, 
£k, arises in a d-dimensional Kondo system, if all conduc- 
tion electron states couple equally to the impurity spin 
It is the length of a finite-size conduction electron 
sea, the "Kondo box", which is so small that its level 
spacing A is comparable to Tk of the bulk system and 



cuts off the logarithmic Kondo correlations. Therefore, 
a box of length £k sustains just one conduction electron 
state within the Kondo scale Tk to form the Kondo sin- 
glet Q, i.e. £k is the size of the Kondo cloud, the Kondo 
screening length. Equating A = Tk, with A the inverse 
of the typical density of states (DOS) in a box of size £k, 
N{tK) = {iK/2TiYSdki-^/{hvp), yields. 



Ik = 27r {UlSdh'^f^) 



l/d 



(1) 



with Sd the surface of the d-dimensional unit sphere 
0. Hence, £k is an intermediate length scale, which 
for d > 2 can be substantially smaller than the coher- 
ence length, 1/kp < £k < £,k, and £k ~ £,k only in 
effectively ID systems. Another length scale, £rkky, 
would arise in dilute Kondo systems as the one when the 
RKKY coupling between neighboring impurities equals 
Tk, tRKKY = [JN{k-p^)Y'HK < Ik 9], where JN{kf) 
is the dimensionless spin coupling. The different physical 
meaning of ^k and £k should be kept in mind for the de- 
sign of related experiments. For example, experiments to 
detect the Kondo cloud via finite system size, like those 
proposed in Refs. , probe £k rather than . These 
experiments should be performed on ID wires in order 
for £k to be in an experimentally accessible range. ID 
Kondo boxes have up to now been realized as ultrashort 
Carbon nanotubes which, however, do not easily per- 
mit persistent or transport current measurements. 

In this work we show numerical evidence that the 
Kondo cloud can be detected via spectroscopy of the 
Kondo resonance in a ID Kondo box. To that end we es- 
tablish a non-trivial proportionality between the Kondo 
spectral weight and the spin screening length as function 
of system size, using large-scale density matrix renormal- 
ization group (DMRG) calculations [l3jll21- The systems 
considered here are ID in the sense that the magnetic im- 
purity is side-coupled to a finite chain of atoms only at 
a single site xq of the chain. This is different from the 
ultrasmall boxes considered in Refs. 0, Q|, where the 
effective hybridization was the same for all states in the 
box. The latter systems may have been realized most re- 



2 



0.04 



"T I ' r 



' I ' I ' I ' I ' I 

128 ^7^N=130 

129- 



-0.1 0.1 0.2 



(B/(D/2) 



-0.1 0.1 0.2 0.3 



0.18 




0.14 



100 



200 







300 



400 



500 



FIG. 1: (Color online) Finite-size and even/odd effects in the 
ID Kondo box. Solid line: Tk as a function of impurity posi- 
tion xo for L — 511, A*' = 512, V = 0.35 and open boundary 
conditions, as obtained from Eq. Open circles: Weight of 
the Kondo peak, Wk, as defined in the text and in the inset. 
The results shown are for even xo only (see text). The inset 
shows the Kondo peak for V — 0.35, a;o = 4 and for various 
successive values of L and A'^. 

cently in molecules 15). As a result of the local coupling 
we observe strong mesoscopic variations of Tk and of the 
spectral features. We analyze, under which mesoscopic 
conditions the above-mentioned proportionality prevails. 

The Hamiltonian for an Anderson impurity with local 
energy Sd and on-site Coulomb repulsion U, side-coupled 
via the hybridization matrix element V to the site xo on 
a ID chain of L sites, reads. 



(2) 



where Hch = -tJ2{^,j),a = 0, ...,L- 1, is 

the free chain Hamiltonian with nearest-neighbor hop- 
ping t > 0. For the evaluations we choose the total elec- 
tron number N near half bandfilling {N = L±2, ej? w 0) 
and use generic parameters for the model in the Kondo 
regime, Sd — —0.55, U — 5, and V as indicated where 
appropriate. All energies are given in units of the half 
bandwidth D = 2t. The Kondo spin coupling is given by 
J ^V^l/\ed\ + l/ied + U)]. 

Tk in finite systems. As mentioned above, for this re- 
alistic model one expects large finite-size effects, because 
the effective impurity-chain coupling, which governs the 
low-energy Kondo physics, depends on the amplitudes of 
the free-electron eigenfunctions of the chain, ^'^(xo), at 
the position xq. The Kondo scale Tk is defined as the 
temperature T at which the 2nd order contribution to 
the spin scattering T-matrix equals the 1st order 0, a 
condition which in the finite system reads, 
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FIG. 2: (Color online) The impurity spectral density Ad{uj) 
for L = 127, iV = 128, V = 0.35 and two different xo; V = 
0.05. Comparison with the non-interacting spectrum {U = 0) 
exhibits the Kondo enhancement of the peak near ef = 0. 
The inset shows the upper Hubbard peak near cu — Sd + U . 

with Ek the levels of the free chain. It is seen that Tk 
itself depends on the impurity position xo |lflllll | and on 
the system size L as well. The strong xq dependence of 
Tk{xq) shown in Fig. ^ is due to the increase of the ID 
local DOS towards the ends of a chain with open bound- 
ary conditions 0|. If xq ist too close to the center of 
the chain (e.g. \xo/L - 1/2| < 160 in Fig. [l|, the log 
contributions in Eq. (|3J are cut off by the level spacing 
of the finite system before the breakdown of perturba- 
tion theory, so that the system stays in the perturbative 
regime for all temperatures, i.e. Tk = (Fig.QJ. Hence, 
in an ultrasmall system the expressions for £k and S.K 
discussed above can, at best, serve to obtain typical val- 
ues for these quantities. We find that the width of the 
Kondo resonance for various Sd, U, V resembles roughly 
Tk of Eq. 0, however obscured by the discreteness of 
the box spectrum. Detecting the Kondo cloud by varying 
the system size then becomes a nontrivial task, since £k 
itself depends on L. Detailed numerical calculations are, 
therefore, needed in order to incorporate these finite size 
effects and to extract the universal features that persist 
under these conditions. 

Numerical method and testing. Applying an efficient 
DMRG code |l^ to the model Eq. lO, we have computed 
the (retarded) impurity Green's function and the equal- 
time spin correlation function at T = 0, 
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respectively, for system sizes of up to L = 511. Here 
Lo is the single-particle excitation energy relative to the 
many-body ground state energy Eq, E — uj + Eq. |0) 
I ^ the DMRG many-body ground state and S'f , the z- 

components of the spin-1/2 operators on the impurity 
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FIG. 3: (Color online) The average C(r) of the spin correla- 
tion function K{r) is shown as function of 1/r for L — 511, 
A'^ = 512 and various hybridization strengths V; r — x — xq, 
xo = 4. Inset: K{r) showing RKKY oscillations. The 
V — 0.55 curve is offset by -0.005 for clarity. 



and on chain site x, r — x ~ xq, respectively. The im- 
purity spectral density is Adai^^) = —-^^i^Gdai^^)- Open 
boundary conditions are applied to facilitate convergence 
of the DMRG algorithm. They also appear appropriate 
for a wire (weakly) coupled to leads. For the dynamical 
quantities we have used both the correction vector (CV) 
method 01 j and the Lanczos method (LM) . For the 
CV method, m = 200 basis states were retained in each 
DMRG iteration, which proved sufRcient to compute the 
residue of the CV (uj + ir]- H)-^di\0) with a precision of 
10"^ for each cu. For the LM we used 3 to 5 target states, 
kept (m — 2600) basis states and carried out 200 Lanczos 
steps to build the Krylov subspace. The comparison of 
the two methods for L up to 128 yields excellent agree- 
ment (better than 0.1 per cent) for lo < Tk and still good 
agreement (better than 10 per cent) even for the highest 
\uj\ K D, where the LM becomes inaccurate. Scaling up 
the system size from L = 128 to i' = 511 reduces the 
frequency range where Lanczos is accurate by a factor 
L/L', which was satisfactory for the calculations in the 
Kondo regime. For the largest systems {L = 511) we, 
therefore, used the numerically less demanding LM. 

Note that all DMRG calculations are done in the 
canonical ensemble with fixed electron number N and 
fixed total spin 5, whereas experimental systems are usu- 
ally coupled to a particle reservoir. Life-time effects of N 
and S are included as a Lorentzian (for the CV method) 
or Gaussian (for the LM) broadening 77 of the energy lev- 
els, with 77 = 0.05 below, xq is chosen near the end of the 
chain, where Tk is large enough (see above and Fig.^ so 
that we can sweep through the crossover from A > Tk to 
A < Tk- Furthermore we choose xq to be even, because 
on all odd sites the chain wave function a,t sp = has a 
node, so that for small L (A > Tk) the impurity would 
be decoupled. 

Results. The T = impurity spectrum Adai'-^) shown 
in Fig. 121 exhibits a rich multiple peak structure even in 
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FIG. 4: (Color online) The weight of the Kondo resonance 
Wk and the numerically determined screening length E*;^ as 
function of size L for xo = 4 and N = L + 1, N odd. 



the single-particle spectral weight near lu ^ e^, induced 
by the discrete local conduction electron spectrum even 
for the largest L, when the impurity is placed close to 
the boundary. The Kondo peak is identified in Fig. |2 as 
the one near uj — through its systematically increas- 
ing weight as the interaction U is switched on, as L is 
increased (Fig. 0}, or as Tk is increased by moving the 
impurity from xq = 50 to = 4 (see also Fig. P). The 
latter would correspond to decreasing T Q in a temper- 
ature dependent measurement. For the local impurity 
coupling V in Eq. (0) we find that the particle number 
parity effect in the position of the spectral features (1 or 
2 peaks within \lu\ < Tk) is essentially washed out by 
finite-size irregularities of the local conduction electron 
spectrum even for small level broadening rj (not shown), 
in contrast to the pronounced even/odd characteristics 
predicted for equal coupling to all conduction states 
However, the even/odd effect is seen in the inset of Fig.Q] 
as an enhancement of the Kondo peak for even as com- 
pared to odd N for fixed system size L. 

The impurity-conduction electron spin correlation 
function K(r), as computed from Eq. (jSJ, is shown in 
the inset of Fig. |3| It displays RKKY oscillations with 
period Xrkky = i^/kp = 2a (a =lattice constant). Its 
overall weight yields s~'^^^ K{r) = — n^, with Ud the 
total impurity occupation number, confirming complete 
screening of the impurity spin, s — 1/2. The average 
C{r) — (K{r) + K{r + l))/2 measures the spin con- 
tent in the Kondo cloud at distance r, while AK{r) = 
\K{r) - K{r + l)|/2 is the amplitude of the RKKY oscil- 
lations. C{r) is shown in Fig.|21 together with the respec- 
tive Ek as calculated from Eqs. (0), The expected 
l/r'' behavior 0] is clearly seen for 1/kp ^ r < £k and 
V = 0.3, 0.34. For smaller £k (V = 0.4, 0.45, 0.55) 
the powerlaw range is too narrow to be observable. For 
r ^ Ik, we find exponential decay, C(r) c>c exp(— 2r/£if) 
(Fig. OJ, and similar for Aif(r). This is expected for the 
correlator of two non-conserved quantities, like Sf , S^, 
with a finite correlation length. In the asymptotic region. 
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FIG. 5: (Color online) Wk like in Fig. El but for A'' even. 

r ^ iK, the exponential behavior should be overidden 
by the slower powerlaw decay, C(r) (x . expected 

from general Fermi liquid arguments [l9l l20l | . The nu- 
merical data show indications of this crossover for the 
largest L and the smallest £k- A more detailed analysis 
of the complex r-dependence will be presented elsewhere. 

For V = 0.3 the £k from Eq. (QJ is 936 > L. C(r) is 
then not cut off by ix but by L, and the conduction elec- 
tron spin density necessary for complete spin screening 
is accumulated at shorter distances, leading to a positive 
y-axis intersection, see Fig. 13 This displays the diffi- 
culty in extracting Ik directly from finite systems and 
the limited applicability of Eqs. iQ, for this purpose. 
Therefore, we combine the results for Adaii-^) and C(r) 
to obtain an experimental signature of the (bulk) screen- 
ing length iji in the finite-size spectra. In doing so one 
must observe that for L < £x, itself becomes size 
and position dependent according to Eqs. (Q, ^ and 
that for our system with fixed total spin there is always 
a total spin 1/2 in the cloud, no matter how small L. 
Therefore, we define the screening length of the finite 
system, £*j^{L), by the volume needed to host a certain 
fraction c of the total spin, J^"^ drK(r) = c, where 
s = 1/2 is the electron spin. The Kondo spectral weight 
is defined as Wk{L) = J^^, dujAdai^^) (cf. Fig. ^ in- 
set), where the boundaries — e', e, are chosen so as to 
cover the Kondo resonance, identified numerically as that 
part of the spectrum around e f which increases as U is 
switched on (cf. Fig. [SJ. The results for both quantities 
are shown in Figs.0I[Slfor c = 0.75 and e = e' = 0.2. For 
odd particle number N (Fig. ^ the nontrivial propor- 
tionahty £*j^{L)/Wk{L) = a{J) for the complete range 
of L is established. We checked that it persists indepen- 
dent of the precise choice of e, e' and c. Both Wk{L) 
and £*x{L) are logarithmically suppressed with decreas- 
ing L (inset of Fig. 0}. For universality reasons we ex- 
pect the proportionality to extend out to L ^ cxd, where 
£*x{L) £k{oo) must saturate at its bulk value. The 
proportionality Wk{L) oc £*k{L) persists for different val- 
ues of J (Fig. 0)), and the corresponding (i) can be 
scaled on top of each other by plotting £*k{L) / £*ii{J) vs 
L/£*j^{J), with the scaling parameter £*j^{J) ~ £k{oo). 
The above relation can be used to determine £*k{L) by 
a spectroscopic measurement and to extrapolate to its 



bulk value, once the proportionality constant a{J) is de- 
termined. Fig. [S] displays for even N , showing an ear- 
lier saturation compared to Fig. 01 as expected from the 
even/odd effect 0. However, we find Wk{L) w const. 
in this case, breaking the above proportionality. By an 
analysis of the spectra this is traced back to the fact 
that for the parameters of Fig. [S] the impurity spectrum 
is dominated by a strong L-independent single-particle 
peak inherited from the free conduction band, while the 
spin structure, C{r), retains its L-dependence. This is 
to emphasize that it is essential to identify the ep peak 
as a Kondo peak first, e.g. by its logarithmic T or i 
dependence, before the above analysis can be applied. 

To conclude, we have analyzed the spectral and the 
spin structure of ultrasmall Kondo systems in the pres- 
ence of strong finite-size fiuctuations and even/odd ef- 
fects using DMRG. Despite these non-universal effects we 
have identified a procedure to measure the spin screen- 
ing length £k by tunneling spectroscopy, e.g. on carbon 
nanotube Kondo boxes. Further research is needed to 
understand the relation £*x{L) — aWxiL) and to deter- 
mine the proportionality factor a(J). 
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